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ABSTRACT 

Brightest cluster galaxies (BCGs) are usually quiescent, but many exhibit star formation. Here we 
exploit the opportunity provided by rest-frame UV imaging of galaxy clusters in the CLASH (Cluster 
Lensing and Supernovae with Hubble) Multi-Cycle Treasury Project to reveal the diversity of UV 
morphologies in BCGs and to compare them with recent simulations of the cool, star-forming gas 
structures produced by precipitation-driven feedback. All of the CLASH BCGs are detected in the 
rest-frame UV (280 nm), regardless of their star-formation activity, because evolved stellar populations 
produce a modest amount of UV light that traces the relatively smooth, symmetric, and centrally 
peaked stellar distribution seen in the near infrared. Ultraviolet morphologies among the BCGs with 
strong UV excesses exhibit distinctive knots, multiple elongated clumps, and extended filaments of 
emission that distinctly differ from the smooth profiles of the UV-quiet BCGs. These structures, 
which are similar to those seen in the few star-forming BCGs observed in the UV at low redshift, are 
suggestive of bi-polar streams of clumpy star formation, but not of spiral arms or large, kpc-scale disks. 

Based on the number of streams and lack of culprit companion galaxies, these streams are unlikely 
to have arisen from multiple collisions with gas-rich galaxies. These star-forming UV structures are 
morphologically similar to the cold-gas structures produced in simulations of precipitation-driven AGN 
feedback in which jets uplift low-entropy gas to greater altitudes, causing it to condense. Unobscured 
star-formation rates estimated from CLASH UV images using the Kennicutt relation range up to 
80 M 0 yr _1 in the most extended and highly structured systems. The circumgalactic gas-entropy 
threshold for star formation in CLASH BCGs at 2 ^ 0.2-0.5 is indistinguishable from that for clusters 
at 2 < 0.2. 

Subject headings: galaxies: clusters: intracluster medium 


1. INTRODUCTION 

The brightest cluster galaxy (BCG) is likely to be the 
dominant galaxy associated with a cluster of galaxies. As 
such, it is often but not always found near the bottom 
of the cluster’s gravitational potential well (e.g., Lauer 
et al. 2014) and experiences a formation and interaction 
history unique to its special position in the hierarchy of 
galaxies (De Lucia & Blaizot 2007). Many BCGs exhibit 
colors and light profiles similar to those of quiescent el¬ 
liptical galaxies, but a significant fraction have extremely 
extended halos, with Sersic indices (Sersic 1968) trending 
closer to n ~ 10 than to the canonical n = 4 expected 
for elliptical galaxies (Caon et al. 1993). 

The unique features and relative rarity of BCGs in typ¬ 
ical galaxy samples have arguably led to some neglect in 
studies of galaxy evolution. However, they provide an 
excellent opportunity to study the relationships between 
star formation, the circumgalactic medium (CGM), and 
feedback from a central active galactic nucleus (AGN). 
Baryons in the CGM are the probable gas supply for 
star formation in galaxies and possibly also for super- 
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massive black-hole growth. Around BCGs the CGM 
is hot enough to study in detail with X-ray observa¬ 
tions, whereas the cooler CGM of less massive galaxies 
is probed only sparsely with UV and optical absorption- 
line studies along lines of sight to bright quasars (e.g., 
Steidel et al. 1994; Tumlinson et al. 2013; Werk et al. 
2014). 

Multi-wavelength observational programs have shown 
that many BCGs are not quiescent systems and exhibit 
activity such as radio emission (> 10 23 W/Hz) (Burns 
1990), extended (^ 10 kpc) emission-line nebulae (Hu 
et al. 1985; Heckman et al. 1989), significant excess blue 
or UV light (Cardiel et al. 1998; Rafferty et al. 2008a; 
Donahue et al. 2010), far-infrared emission from warm 
dust and PAHs (Quillen et al. 2008; Donahue et al. 2011), 
and vibrationally excited molecular hydrogen (Donahue 
et al. 2000). This enhanced activity in the BCG is 
strongly correlated with the thermodynamic state of the 
hot CGM surrounding the BCG (e.g., Cavagnolo et al. 
2008; Rafferty et al. 2008a; Cavagnolo et al. 2009; Sun 
2009). In rich clusters of galaxies, all such active BCGs 
are embedded in cores or coronae of dense, low-entropy, 
high-pressure gas, commonly called “cool cores” and for¬ 
merly known as “cooling flows” (Donahue & Voit 2004). 
The radiative cooling time of the hot gas surrounding 
active BCGs is typically < 1 Gyr. 

While 10-30% of BCGs in optically-selected samples 
exhibit either radio-loud AGN or star-formation activ¬ 
ity or both (e.g., von der Linden et al. 2007; Edwards 
et al. 2007), the fraction rises to > 70% in the local sam- 
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pie of X-ray selected clusters (Crawford et al. 1999) and 
in a representative sample of X-ray clusters (Donahue 
et al. 2010). Some of these BCGs have among the high¬ 
est SFRs in the £ < 0.6 universe (e.g., McNamara et al. 
2006; McDonald et al. 2013). The relationships observed 
between molecular gas content and star-formation rate in 
star-forming BCGs are remarkably similar to those ob¬ 
served in star-forming galaxies and starbursts. Molecular 
gas “depletion times” based on molecular hydrogen (H 2 ) 
masses derived from CO line luminosities (Edge 2001) 
and star-formation rates derived from FIR emission are 
^ 10 9 years for both star-forming galaxies and BCGs 
with SFRs of a few solar masses per year, and drop to 
~ 10 7 — 10 8 years for both starburst galaxies and BCGs 
with SFRs of tens or more solar masses per year (Voit & 
Donahue 2011; O’Dea et al. 2008). 

One hypothesis for the link between a low-entropy 
CGM and activity in the BCG is that the lowest en¬ 
tropy gas in the CGM cools and provides cold gas to fuel 
both star formation and the AGN. A long-standing mys¬ 
tery is that unmoderated radiative cooling of the ICM 
would produce star formation rates at least an order 
of magnitude higher than observed (e.g., Fabian et al. 
1991), hence the “cooling flow” problem. Spectroscopic 
evidence from the XMM-Newton grating spectrometer 
conclusively demonstrated that the ICM is not radia- 
tively cooling from X-ray emitting temperatures (Peter¬ 
son et al. 2003), leading most to assume that some kind 
of heating process involving the AGN nearly balances 
radiative cooling (McNamara & Nulsen 2007, 2012). A 
second mystery is that, despite the strong correlation 
between low-entropy hot gas and cold, molecular, star¬ 
forming gas, it is not at all clear how hot, dust-free, X-ray 
emitting gas becomes dusty, cold molecular gas (Voit & 
Donahue 2011). 

Observations of BCGs and the surrounding CGM sug¬ 
gest that the relationships between the AGN, its kinetic 
energy output, the BCG’s star-formation activity, and 
the thermodynamic state of the intergalactic gas must 
be fairly cozy. Simulations of the broader population of 
galaxies indicate that star formation in the most massive 
galaxies must somehow be quenched by AGN feedback, 
because without it massive galaxies are not only much 
more luminous than galaxies we see today, but bluer as 
well, owing to high levels of continuous star formation 
(e.g., Saro et al. 2006). This problem is rectified in sim¬ 
ulations that incorporate AGN feedback, indicating that 
AGNs are the likely culprits quenching continuous star 
formation in massive galaxies. The phenomena we ob¬ 
serve in cool-core BCGs are likely to be signposts of AGN 
feedback in action. These AGN do not eliminate star for¬ 
mation, but are energetically capable of preventing the 
extreme star formation that would ensue if the ICM were 
allowed to cool unabated into molecular clouds (e.g., Fab- 
jan et al. 2010; Dubois et al. 2013; Ragone-Figueroa et al. 
2013; Martizzi et al. 2014). 

We and others have recently proposed a framework 
for precipitation-driven AGN feedback that account for 
these tight relationships (e.g., McCourt et al. 2012; 
Sharma et al. 2012; Gaspari et al. 2012; Voit & Donahue 
2015; Voit et al. 2015). In this framework, radiative cool¬ 
ing of thermally unstable gas causes condensation and 
precipitation of cool clouds out of the hot CGM, which 
then feed star formation and fuel growth of the central 


supermassive black hole. The criterion for precipitation, 
indicated by numerical simulations of feedback and radia¬ 
tive cooling in a spherical geometry is for the cooling time 
to be shorter than about ten times the free-fall time. In 
the highest mass galaxies, precipitation causes a jump in 
jet power from the AGN, which subsequently limits cool¬ 
ing. The jets are also important for lifting low-entropy 
gas from the BCG’s center out to larger radii where it 
becomes unstable to condensation and creates filaments 
and knots of cooler gas. If dust grains released by the 
BCG’s old stellar population can survive until they are 
uplifted and incorporated into the condensing gas, then 
they can help to promote cooling and can also serve as 
nucleation sites for converting a greater proportion of re¬ 
fractory elements into solid form. Numerical simulations 
have begun to reproduce observable features of this phys¬ 
ical model (e.g. Li & Bryan 2014a, b). 

Here we present UV images from the Hubble Space 
Telescope Multi-Cycle Treasury program CLASH (Post¬ 
man et al. 2012b) that illustrate the most recent loca¬ 
tions of star formation in active BCGs. Very few high- 
resolution UV images exist for brightest cluster galax¬ 
ies, and the handful available at z < 0.3 have revealed 
filamentary, clumpy structures (Koekemoer et al. 1999; 
O’Dea et al. 2004, 2010; Hicks et al. 2010; Oonk et al. 
2011). CLASH has substantially boosted the scant high- 
resolution UV coverage of BCGs by collecting high- 
quality images of 25 clusters of galaxies at 0.2 < z < 0.9 
through 16 filters spanning the ultraviolet to the near- 
infrared. The rest-frame UV images of these galaxy clus¬ 
ters have no precedent in sensitivity or spatial resolution 
at similar redshifts. We briefly describe the CLASH clus¬ 
ter sample in §2. We describe the data analysis of the UV 
images for the 25 CLASH BCGs, estimate unobscured 
UV star-formation rates, and analyze UV morphology 
in §3. We discuss these results and compare the UV 
morphologies to simulations of star formation in central 
cluster galaxies experiencing AGN-jet feedback in §4. We 
summarize our results in §5. 

Throughout this paper we assume cosmological param¬ 
eters of Dm = 0.3, Da = 0.7, and Hq = 70 I 170 km s -1 
Mpc -1 . 

2. CLASH 

The Cluster Lensing and Supernova with Hubble 
(CLASH) Multi-Cycle Treasure (MCT) program (Post¬ 
man et al. 2012b) utilized 525 orbits (over ^ 500 hours 
of observing time) from 2011-2013 (HST Proposal ID 
12456). The sample includes 25 clusters of galaxies with 
z = 0.2 —0.9 and global X-ray temperatures greater than 
5-6 keV. Twenty of these clusters were chosen to be ex¬ 
amples of relaxed, massive systems, characterized by reg¬ 
ular X-ray morphology on radial scales greater than a 
few hundred kpc, with BCGs that are well-centered and 
aligned with the X-ray centroids. A significant number 
of these clusters were included in the Mantz et al. (2008) 
study of relaxed X-ray systems. The remaining five clus¬ 
ters were added to the sample because of their lensing 
properties and ability to magnify the high-redshift uni¬ 
verse (Zheng et al. 2012; Coe et al. 2013). All of the 
CLASH clusters have Chandra observations, and 15 were 
also observed with XMM (Donahue et al. 2014). 

One of the unprecedented features of the CLASH study 
was its use of 16 broadband filters, spanning UV to near- 
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IR wavelengths. The motivation for this broad coverage 
was to significantly improve the photometric redshift es¬ 
timates of lensed galaxies, but a side benefit is that the 
cores of these clusters received UV exposure otherwise 
rather difficult to obtain. Only a few of the CLASH 
BCGs had published UV observations prior to CLASH. 
A383, mlll5, r2129, and ms2137 have GALEX-based 
SFR upper limits of < 2 — 3 Mq yr -1 . Only one, ml720, 
had a published estimate of its UV star-formation rate 
(2.7 ± 0.7 M© yr- 1 ; Hoffer et al. 2012). 

These UV images provide our first high resolution view 
of a sample of BCGs in massive clusters at moderate 
redshifts, but many of them have been the subject of 
studies at other wavelengths, only a few of which we 
have space to mention here. Several of the BCGs are 
known to be mid-IR sources of dust emission, some of 
which can be produced by radiative heating of dust by 
light from evolved stars, and some show clear evidence 
for recent star formation. Four CLASH BCGs have mid- 
IR detections with Spitzer (a209, A383, ms2137, r2248 
and rl532), and there are mid-IR upper limits for several 
more (cl226, A2261, and r2129). The most extreme star¬ 
forming system known in CLASH prior to this study was 
rl532, with a Spitzer-based SFR of 110 M© yr -1 (Hoffer 
et al. 2012). The most extreme AGN-infhienced BCG is 
arguably ml931 (Ehlert et al. 2011). Uniformly derived 
profiles of temperature, ICM electron density, gas en¬ 
tropy, gas mass, and total mass based on analysis of the 
pressure profiles are provided in Donahue et al. (2014) 
for all 25 of the CLASH clusters. Gravitational lensing 
profiles for a significant fraction of the CLASH sample 
are published (Umetsu et al. 2014; Merten et al. 2014; 
Zitrin et al. 2014). All of these sources have either up¬ 
per limits or detections from radio surveys, which we will 
discuss in § 4.3. 

3. CLASH OBSERVATIONS 

The CLASH clusters, mean redshifts, observation 
dates, and filter names are given in Table 1. Filter 
names are provided using the convention that an HST 
filter named FnnnW is listed as “nnn” in the table. For 
HST, the filters shortward of 1 micron are named for the 
approximate central wavelength of the filter in nanome¬ 
ters; filters longward of 1 micron (which we call infrared 
or IR for short) are named for this wavelength in units of 
10 nanometers (so 1.6 micron would correspond to 160 
on this scale). 

3.1. HST Data Reduction and Photometry 

The HST ACS and WFC3 imaging data for these clus¬ 
ters were processed through initial detector-level calibra¬ 
tion to remove instrumental signatures from each expo¬ 
sure. Exposures in all filters were subsequently aligned 
onto a uniform astrometric grid to an accuracy better 
than ~ 2-3 milliarcseconds, passed through cosmic-ray 
and bad-pixel rejection, and combined into final mosaics 
for each filter at pixel scales of 30 milliarcseconds and 
65 milliarcseconds. All of this was done using a modified 
version of the MosaicDrizzle pipeline (Koekemoer et al. 
2011, 2002). Further details are presented in the CLASH 
overview paper (Postman et al. 2012b). We used mosaics 
binned to 65 milliarcseconds per pixel for this study. For 
the WFC3 UV images, later versions of our mosaics were 
significantly improved over the first version of publicly 


available mosaics because of better flat-fielding and sky 
treatment. For the sake of reproducibility, we provide 
the mosaic production date in Table 1. 

Rest-frame UV was defined to be 280 nm in the rest- 
frame of the BCG. At this wavelength, the UV from re¬ 
cently formed hot stars, if present, readily exceeds the 
light emitted by older stars. At shorter wavelengths, the 
relative contribution of “UV-upturn” stars increases be¬ 
low 250 nm (e.g., Brown 2004) with considerable vari¬ 
ation from one galaxy to the next, even in a sample of 
quiesent ellipticals (e.g., Rich et al. 2005; At lee et al. 
2009). To allow an estimate of the level of the evolved 
stellar population, we chose the image obtained at a rest 
wavelength of ^ 1 micron, chosen because 1 micron is 
approximately the peak of the stellar-light spectrum for 
passive galaxies and evolved stars. For the highest red- 
shift CLASH clusters, cl226 and m0744, that wavelength 
was not in the CLASH coverage, so the near-IR (NIR) 
photometry is based on the reddest bandpass available 
(F160W). Resolved dust maps, emission-line maps, and 
optical long-slit spectroscopy will be presented in a future 
work (Fogarty et al., in prep.) For accurate color pho¬ 
tometry, the UV images were degraded with a Gaussian 
kernel to match the PSF in the NIR prior to photometry 
measurements. 

We measured the total UV emission from each BCG 
using the IRAF 2.15.1a package apphot within a circular 
aperture centered on the BCG 6 , with coordinates and 
radii listed in Table 2. For BCGs with structures visi¬ 
ble in the UV, the aperture size was chosen to provide 
an estimate of the total amount of UV light from the 
highest surface brightness structures in the BCG. If the 
field was uncrowded and the UV source centrally concen¬ 
trated and extended, a default aperture of 14 kpc (equal 
to 10 /h kpc) was chosen to estimate the color of quies¬ 
cent BCGs. However, this aperture size was reduced in 
crowded fields to avoid including light from neighboring 
galaxies, or in the case of ml 149, a prominent and ex¬ 
tended lensed star-forming background galaxy. UV and 
1.0-micron photometry are based on identical apertures 
for the same galaxy. 

To investigate the effects of large-scale scattered light, 
we obtained UV background measurements in annuli of 
multiple sizes surrounding the aperture. Robust UV pho¬ 
tometry measurements are optimized with an inner ra¬ 
dius for the background annulus of about 15", although 
the result is fairly insensitive to choices smaller and larger 
than this. We further investigated background systemat- 
ics by randomly moving the aperture/annulus combina¬ 
tion for each cluster, avoiding sources of obvious emission 
or scattered light. We then performed aperture photome¬ 
try on these blank pieces of sky in the same manner as on 
the BCG. From these measurements, we found the stan¬ 
dard deviation in flux, using a 3-cr cut to exclude possibly 
real but low-flux UV sources. We use this standard de¬ 
viation as the systematic uncertainty in our photometric 
measurements resulting from scattered light in the field 
and other flat fielding issues. 

The NIR sky background was estimated in an aper¬ 
ture well outside the brightest part of the BCG to give 
our best estimate of the ambient sky. The background 
sky levels were based on a median estimate with 3-cr 

6 http://iraf.noao.edu 


4 


Donahue et al. 


Table 1 

HST Observations and Processing Dates 


Full Name 

Cluster 

z 

ObsDate (UV) 

ProcDate 

UV 1 

texp 

(?) 

NIR 1 

texp 

(s) 

ObsDate (NIR) 

Auv 

mag 

Anir 

mag 

Abell 1423 

A1423 

0.213 

12/28/12 

2/27/13 

336 

4890 

125 

2514 

12/31/12 

0.085 

0.014 

Abell 209 

A209 

0.206 

9/9/12 

1/19/13 

336 

4752 

125 

2514 

8/12/12 

0.086 

0.014 

Abell 2261 

A2261 

0.224 

4/3/11 

1/23/13 

336 

4817 

125 

2514 

3/9/11 

0.193 

0.031 

Abell 383 

A383 

0.187 

12/28/10 

7/2/11 

336 

4868 

125 

3320 

1/5/11 

0.14 

0.023 

Abell 611 

A611 

0.288 

3/1/12 

6/18/12 

336 

4782 

125 

2514 

1/28/12 

0.255 

0.042 

WARP 1226.9+3332 

cl226 

0.890 

5/24/13 

6/26/13 

475 

4396 

160 

5129 

5/9/13 

0.063 

0.014 

MACSJ0329.6-0211 

m0329 

0.450 

9/26/11 

11/5/11 

390 

4902 

160 

5029 

8/18/11 

0.235 

0.031 

MACS J0416.1-2403 

m0416 

0.395 

8/18/12 

1/19/13 

390 

4814 

160 

5029 

7/24/12 

0.16 

0.021 

MACSJ0429.6-0253 

m0429 

0.399 

12/11/12 

2/27/13 

390 

4562 

140 

2312 

12/5/12 

0.235 

0.037 

MACS J0647.8+7015 

m0647 

0.584 

11/16/11 

12/2/11 

435 

4248 

160 

5229 

10/22/11 

0.402 

0.057 

MACSJ0717+3745 

m0717 

0.548 

10/29/11 

12/30/11 

435 

4052 

160 

5029 

10/10/11 

0.277 

0.039 

MACSJ0744.9+3927 

m0744 

0.686 

11/17/11 

12/30/11 

475 

4022 

160 

5029 

10/29/11 

0.189 

0.03 

MACSJ1115.8+0129 

mlll5 

0.355 

1/31/12 

3/5/12 

390 

4554 

140 

2012 

1/7/12 

0.151 

0.024 

MACSJ1149.6+2223 

mll49 

0.544 

2/13/11 

3/14/11 

435 

3976 

160 

5029 

1/16/11 

0.083 

0.012 

MACS J 1206.2-0847 

ml206 

0.439 

5/25/11 

8/15/11 

390 

4959 

160 

5029 

1/16/11 

0.246 

0.032 

MACSJ1311.0-0311 

ml311 

0.494 

6/12/13 

7/9/13 

435 

4172 

160 

2011 

4/22/13 

0.112 

0.016 

MACSJ1423.8+2404 

ml423 

0.545 

2/5/13 

3/13/13 

435 

4196 

160 

5029 

1/19/13 

0.099 

0.014 

MACSJ1720.2+3536 

ml720 

0.387 

5/22/12 

6/28/12 

390 

4810 

140 

2312 

5/9/12 

0.147 

0.023 

MACSJ1931.8-2635 

ml931 

0.352 

5/28/12 

6/28/12 

390 

4620 

140 

2312 

4/14/12 

0.431 

0.068 

MACSJ2129-0741 

m2129 

0.570 

6/24/11 

1/22/12 

435 

3728 

160 

5029 

5/15/11 

0.273 

0.039 

MS 2137.3-2353 

ms2137 

0.313 

9/29/11 

1/21/13 

390 

4827 

140 

2312 

9/9/11 

0.198 

0.031 

RXJ1347.5-1145 

rl347 

0.451 

6/22/11 

8/12/11 

390 

4820 

160 

7741 

5/21/10 

0.241 

0.032 

MACSJ1532.8+3021 

rl532 

0.363 2 

3/29/12 

4/16/12 

390 

4840 

140 

2312 

3/16/12 

0.115 

0.018 

RXJ2129.6+0005 

r2129 

0.234 

5/31/12 

10/16/12 

336 

4574 

125 

3420 

4/3/12 

0.18 

0.029 

ABELL S1063 

r2248 

0.348 

9/24/12 

1/19/13 

390 

4740 

140 

2312 

9/12/12 

0.048 

0.007 


1 Filter names are listed here in a compressed format, such that nnn represents filter FnnnW for the ACS or WF3. 

2 Redshift of this cluster and its BCG is based on the SDSS emission-line spectrum of the BCG. 


clipping. After aperture fluxes were converted to AB 
magnitudes using the ACS and WFC3 zeropoints, 7 all 
photometry was corrected for Galactic extinction at the 
observed bandpass (Table 1), ranging from 0.08 — 0.45 
mag in the UV and < 0.01 — 0.07 mag in the NIR, us¬ 
ing HST bandpass-specific extinction corrections assum¬ 
ing Ry =3.1 and based on (Fitzpatrick 1999)based on 
Schlafly & Finkbeiner (2011), which is a recalibration of 
Schlegel et al. (1998), and. 

We estimated K-corrections based on the assumption 
that the NIR spectrum is dominated by a 10 Gyr old 
stellar population, which was modeled with HST band- 
passes convolved with Starburst99 models (Leitherer 
et al. 1999). The estimated K-correction is less than 0.1 
magnitudes for all the clusters and does not affect our 
results in an interesting way. We did not attempt to ho¬ 
mogenize the UV photometry for small offsets from the 
nominal rest-frame wavelength of 280 nm, since the UV 
light from a population of hot stars is nearly flat in F v . 
We verified this assumption by comparing the 280 nm 
rest-frame photometry with bluer data, when available. 

Statistical errors for both UV and NIR were computed 
by apphot based on counting statistics, obtained by set¬ 
ting the appropriate electron-to-count-rate ratio equal 
to the effective gain, and the standard deviation in the 
background annuli. (The mosaics in the CLASH pipeline 
have the units of electrons sec -1 , so to compute the ap¬ 
propriate counting statistics for a given mosaic, the gain 
parameter in the apphot task was set to the exposure 
time.) Because the galaxies are so bright in the NIR in 
the aperture of interest, the exact NIR background and 
its uncertainty are negligible for this work. 

7 http://www.stsci.edu/hst/acs/analysis/zeropoints/zpt. 
py , http://www.stsci.edu/hst/wfc3/phot_zp_lbn 


Systematic and statistical flux errors were added in 
quadrature to estimate the la uncertainties given in Ta¬ 
ble 2. Color uncertainties include an estimated calibra¬ 
tion uncertainty of 2%, to account for uncertainties in 
absolute calibration, relative calibration between the two 
filters, and the unknown underlying spectral shape com¬ 
pared to that assumed to estimate the AB magnitude at 
a given wavelength. 

3.2. UV Star-Formation Rates 

In order to estimate the UV light from recent star for¬ 
mation, we first need to account for the UV light con¬ 
tributed by the evolved stellar population. For this pur¬ 
pose, we use the quiescent BCGs as templates. The fact 
that the BCGs are at different redshifts introduces some 
scatter into our estimate, but not enough to affect these 
estimates significantly. This method has the advantage of 
obtaining the UV contribution from evolved stars from 
an otherwise similar population of galaxies: the BCG 
population itself. 

The UV-NIR colors of quiescent BCGs are fairly well 
behaved within the CLASH sample. We find an average 
280-1000 micron color of 5.13 ±0.35 (omitting ml206 be¬ 
cause of contamination by a lensed galaxy). This is con¬ 
sistent with the mean and dispersion of quiescent BCGs 
in other studies, and a baseline (established by the 3 
reddest BCGs) of 5.5 ± 0.5. In Hoffer et al. (2012) we 
found that the typical UV-K color of relatively nearby, 
quiescent BCGs was 6.6 ± 0.3, and typical J — K colors 
are ^ 1. The BCGs with the faintest inferred UV lumi¬ 
nosities also have the reddest UV-1.0 micron colors of 
^5 — 6, consistent with UV-NIR estimates from Hoffer 
et al. (2012); Hicks et al. (2010); Donahue et al. (2010). 
We use the average color to correct for the UV produced 
by the evolved stars in the BCG. Using a redder color 
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Table 2 

UV and NIR Photometry 


Cluster 

Radius 

Radius 

UV (AB) 

± 

UV-IR 

± 

Log L uv 

SFR 

SFR 

SFR 

K 0 C 


arcsec 

kpc 

mag 

mag 

mag 

mag 

erg s —1 Hz -1 

M G yr” 1 

corr 

err 

keV cm 2 

ml931 

5.75 

28.5 

18.52 

0.03 

2.04 

0.04 

29.77 

83.08 

78.12 

2.34 

14 ±4 

rl532 

7.02 

35.5 

19.18 

0.03 

2.83 

0.04 

29.54 

48.62 

42.61 

2.64 

17 ±2 

m0329 

3.25 

18.7 

20.22 

0.02 

3.3 

0.03 

29.35 

30.99 

25.08 

2.39 

11 ±3 

m0429 

3.90 

20.9 

20.38 

0.04 

3.75 

0.05 

29.16 

20.11 

14.31 

2.07 

17 ±4 

ml423 

1.20 

7.7 

21.39 

0.01 

3.14 

0.02 

29.08 

16.70 

13.95 

1.15 

10 ±5 

rl347 

2.50 

14.4 

20.91 

0.02 

3.81 

0.03 

29.07 

16.50 

11.47 

1.75 

12 ±20 

mlll5 

1.53 

7.6 

21.33 

0.01 

3.38 

0.02 

28.66 

6.37 

5.06 

0.52 

15 ±3 

ms2137 

1.56 

7.2 

21.16 

0.02 

4.07 

0.03 

28.60 

5.56 

3.41 

0.66 

15 ±2 

ml720 

2.70 

14.2 

21.6 

0.05 

4.54 

0.05 

28.64 

6.08 

2.45 

0.74 

24 ±3 

A383 

2.30 

7.2 

20.46 

0.03 

4.36 

0.04 

28.37 

3.29 

1.63 

0.41 

13 ±2 

r2129 

3.80 

14.2 

21.16 

0.09 

4.98 

0.09 

28.31 

2.86 


< 0.4 

21 ±4 

m0744 

2.00 

14.2 

22.73 

0.13 

4.6 

0.13 

28.78 

8.50 


< 3.1 

42 ± 11 

ml311 

2.41 

14.6 

22.29 

0.14 

4.69 

0.14 

28.61 

5.75 


< 1.9 

47 ±4 

A2261 

3.93 

14.2 

20.9 

0.06 

5.47 

0.07 

28.37 

3.29 


< 2.8 

61 ±8 

ml206 

2.50 

14.2 

21.81 

0.04 

4.47 

0.05 

28.68 

6.75 


< 3 a 

69 ± 10 

A1423 

4.10 

14.2 

21.22 

0.13 

4.96 

0.13 

28.19 

2.19 


< 0.4 

68 ± 13 

a209 

2.50 

8.5 

21.79 

0.11 

5.48 

0.11 

27.93 

1.20 


< 1.1 

106 ± 27 

A611 

1.69 

7.3 

22.93 

0.13 

5.69 

0.14 

27.81 

0.90 


< 1.7 

125 ±18 

cl226 

0.78 

6.1 

24.68 

0.17 

5.37 

0.17 

28.28 

2.67 


< 1.5 

166 ± 45 

r2248 

1.24 

6.1 

22.39 

0.04 

4.91 

0.04 

28.21 

2.29 


< 0.5 

170 ± 20 

m0416 

2.60 

13.9 

22.26 

0.09 

5.24 

0.10 

28.39 

3.48 


< 0.8 

400 ± 100 

m0647 

1.11 

7.3 

23.8 

0.29 

5.19 

0.29 

28.19 

2.14 


< 0.3 

225 ± 50 

m0717 b 

1.95 

12.5 

22.63 

0.15 

4.85 

0.15 

28.59 

5.40 


< 1.4 

220 ± 96 

mll49 

1.11 

7.0 

23.63 

0.15 

4.73 

0.15 

28.18 

2.11 


< 0.7 a 

280 ± 40 

m2129 

1.10 

7.2 

24.03 

0.22 

5.07 

0.22 

28.07 

1.64 


< 0.1 

200 ± 100 


Note. — The CLASH clusters with Kq < 30 keV cm 2 are in the top section, ordered by estimated UV star formation 
rate, from large to small. The next section are clusters with Kq > 30 keV cm 2 , ordered by estimated core entropy Kq. The 
bottom section lists the five high magnification clusters from CLASH, in order of right ascension. All of the information 
in this table is from this work except for the Kq estimates, which are primarily from ACCEPT (Cavagnolo et al. 2009), 
except for m0416 and m2129 (Donahue et al. 2014). 

a The clusters ml206 and ml 149 both have known lensed background systems that contaminate the UV images of their 
BCGs. The upper limit on excess UV in ml206 associated with star formation or other BCG activity is set equivalent to 
the level of detected excess, while the aperture used for ml 149 avoids the lensed galaxy. 

b The cluster m0717 does not have a prominent BCG or even a double BCG. We have identified a bright candidate BCG 
with a spectroscopic redshift consistent with being a cluster member, as discussed in the text. Coordinates for all BCGs, 
along with GALEX photometry for comparison, are provided in Table 4 

c The estimated central entropies Kq = Tx/n in units of keV cm 2 are from ACCEPT (Cavagnolo et al. 2009). For 
the non-ACCEPT cluster m0416, the Kq is based on electron density and temperature profiles from the X-ray analysis 
described in Donahue et al. (2014). M0416 is an interacting cluster, so the high gas entropy in the center is not surprising. 


would increase our SFR estimates most significantly for 
the BCGs with the lowest rates. 

We convert excess UV luminosity to an approximate 
star-formation rate using the UV conversion from Ken- 
nicutt (1998). We note that the Kennicutt SFR conver¬ 
sion is based on the Salpeter (1955) initial mass function 
(IMF) and the assumption of a nearly constant star for¬ 
mation rate over the last 100 million years or so. There 
may be differences in the low-mass end of the IMF be¬ 
tween BCGs and the star-forming galaxies studied by 
Kennicutt (see e.g., van Dokkum & Conroy 2010), but 
for our purposes her the relative SF rates are more rele¬ 
vant than the absolute rates. For the BCGs which may 
be relatively young starbursts, these estimated UV rates 
underestimate the unabsorbed SFR. We have not at¬ 
tempted to estimate the star formation that is obscured 
by dust clouds, but typically, mid-IR based SFR are sim¬ 
ilar to or larger than the unobscured SFRs (e.g., Hoffer 
et al. 2012). 

We find that of the 25 CLASH clusters, seven ex¬ 
hibit unambiguous UV excesses, corresponding to unab¬ 
sorbed star-formation rates of 5-80 solar masses per year 
(ml931, rl532, m0329, m0429, ml423, rl347, mlll5). If 
we adopt a detection threshold of 3 sigma, three other 


clusters exhibit UV excesses possibly associated with star 
formation (A383, ml720, and ms2137). Of these three, 
only ms2137 fails to show morphological evidence for a 
young stellar population clearly distinct from the evolved 
stellar population and could contain a low-luminosity 
AGN. 

Four of the CLASH BCGs have Herschel mid-IR detec¬ 
tions from Rawle et al. (2012), from which they derive 
SFRs of 46 ± 3, 4 ± 0.2, 1.6 ± 0.1, and 1.7 ± 0.1 M e 
yr -1 for ml423, A383, rl720, and r2129, respectively. 
These quantities are quite similar to our UV estimates 
and limits. We did not detect excess UV emission from 
r2129, but distinguishing between excess UV from young 
stars and the UV emission from an evolved population 
becomes difficult with photometry alone at ^ 1 M© yr -1 . 
Ten other CLASH BCGs in Rawle et al. (2012) were not 
detected in the mid-IR. Among these IR non-detections, 
we find a UV excess for only ms2137. 

In the Hoffer et al. (2012) compilation of SFRs, a209 
has an estimated SFR of 0.9 ±0.1 M 0 yr -1 , based on 
24 micron photometry, while our measurements give an 
upper limit of < 1.1 M© yr -1 . The Hoffer SFR for A383 
of 1.6 ± 0.2 is similar to our UV rate of 1.6 ± 0.4 M© 
yr -1 . The UV-quiescent BCG in r2248 (also known as 
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AS1036) has a small IR excess corresponding to 0.8T0.13 
M 0 yr _1 . As mentioned above, rl532 is a known star¬ 
bursting BCG with a mid-IR SFR of 110 ± 22 M 0 yr -1 
and a UV rate of ^ 43 M 0 yr -1 . In general, mid-IR 
observations around 100 microns, near the peak of the 
dust emission, are more sensitive to low levels of star 
formation or other activity because of the uncertain UV 
contribution from evolved stars. However, dust heated 
by light from the old stellar population can also produce 
FIR emission, so low-level SFRs for BCGs based on sin¬ 
gle IR photometry points can also be inaccurate (Don¬ 
ahue et al. 2011). Ideally, one would like more accurate 
measures of star formation based on multi wavelength ob¬ 
servations. However, in most of the star-forming CLASH 
BCGs, we are not hampered by uncertainties in UV flux 
from the evolved population because those galaxies ex¬ 
hibit both UV excesses of clear statistical significance 
and morphological indications that the bulk of the UV 
emission is not coming from the evolved stellar popula¬ 
tion. 


3.3. UV Morphologies 

We show raw (un-smoothed) UV images of the star¬ 
forming CLASH BCGs in Figure 1, and smoothed UV 
images of the rest of the CLASH BCGs in Figure 2. 
Two of the clusters with quiescent BCGs, ml206 and 
ml 149, have well-known lensed background features that 
contaminate UV estimates from the core of the cluster. 
m0744 and ml311 are somewhat bluer than the average 
quiescent BCG, but they both show smooth isophotes 
and little morphological evidence for a distinct star¬ 
forming population. 

In order to characterize the UV and IR morphologies 
of the BCGs in the 20 X-ray selected CLASH clusters we 
compute low-order moments of their surface brightness 
distributions. For each BCG we compute these moments 
inside a circular region that excludes other galaxies and 
estimate a scale and position angle from the resulting 
moment matrix. This exercise can be done analytically, 
so we provide the full procedure below for completeness. 

The surface brightness centroid (C x , C y ) is deter¬ 
mined by the normalized sums C x = ^i(x,y)x and 
C y = ^i(x,y)y, where I(x,y) is the surface brightness 
of pixel (x, y) and i(x, y) = /(x, y)/^2l(x, y). The sums 
are over all (x, y) inside a circular aperture with a ra¬ 
dius chosen to avoid excess off-axis contamination from 
other galaxies. Small variations in this aperture size 
do not change the result. The matrix elements of the 
next moment of the surface brightness distribution are 
M X x = i(x, y) * {x C x ) , Myy = ^ i(x, y) * (y Cy) 
and M xy = M yx = ^i(x,y) * (x - C x )(y - C y ). These 
terms define a symmetric 2x2 matrix 

*-(u) <» 


with M xx = a, Myy = 6, and M xy = M yx = c. This 
matrix can be diagonalized in a standard way with a 
rotation corresponding to 


tan 77 = 


a — b 


+ 1 


1/2 


a — b 
2c 


( 2 ) 


giving the eigenvalues 


A± 


a + b ^ 
2 


a — b 
2c 


2 



( 3 ) 


as diagonal elements. These elements represent one way 
to quantitatively estimate the extent of a complex region. 
We define the terms x s = and y s = y // W. The 

position angle 77 (PA in Table 3) is defined to be that 
of the semi-major axis, in degrees east of north (where 
north is aligned along the y axis as a convention in our 
mosaics), and axis ratio is defined to be the minimum of 
y s /x s or x s /y 3 . 

We computed these quantities for both the UV and 
NIR images and report the results in Table 3. The clus¬ 
ters that host BCGs with clear evidence for star forma¬ 
tion are listed first in this table, approximately in order 
of star-forming region color, from blue to red. All of the 
surface brightness moments for the star-forming BCGs 
were computed within a region with a radius of about 
5", while the moments computed for quiescent BCGs 
were restricted to regions of 1.5 — 3.5" to avoid includ¬ 
ing an unnecessarily large number of sky- or readnoise- 
dominated pixels and to avoid confusion with unrelated 
UV sources. However, the UV and NIR regions used 
to compute morphological quantities for each cluster are 
exactly the same, so that comparison between results 
from the UV and the NIR from cluster to cluster is rel¬ 
evant. The uncertainties for all quantities in this ta¬ 
ble were computed by calculating the standard deviation 
of morphological parameters using 100 bootstrapped im¬ 
ages with gaussian noise proportional to the inverse of 
the square root of the mosaic weight matrix. The sta¬ 
tistical uncertainties for properties derived from the IR 
image were negligible compared to those derived for the 
UV images. 


4. DISCUSSION 
4.1. Star Formation Rates 

Table 2 gives unobscured star formation rates based 
on UV photometry for CLASH BCGs. We have made an 
empirical estimate of the amount of UV associated with 
the evolved stellar population, using the red UV-IR col¬ 
ors of the quiescent BCGs in our own sample. These col¬ 
ors are consistent with those found in previous work on 
low-redshift UV emission from BCGs (Hicks et al. 2010; 
Donahue et al. 2010; Hoffer et al. 2012). Our reported 
UV star formation rates have been adjusted for this con¬ 
tribution assuming a mean UV-NIR color for BCGs of 
5.5. The uncertainty in this component translates to a 
minimum systematic uncertainty in the star formation 
rate of ~ 0.3 M 0 yr -1 for these galaxies. 

As we noted previously, the UV luminosity is a direct 
measurement of the unobscured star formation rate. But 
of course much of the star formation in a galaxy can be 
obscured. Galaxy studies indicate that for typical star¬ 
forming galaxies such as the Milky Way, the unobscured 
star formation rate is approximately equal to the ob¬ 
scured star formation rate, such that SFR^/y ~ SFR fir- 
Typically, the unobscured (UV) star formation rate for a 
galaxy corresponds to 10-25% of the total amount of star 
formation in galaxies with SFR > 2M 0 yr _1 , based on 
comparisons with SFRs estimated from dust emission in 
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Figure 1 . The top 9 images show the central 50 kpcx50 kpc of the 280 nm broad-band HST images for star-forming CLASH BCGs. A 
blue bar on each frame indicates a angular scale. Images are oriented with North up and East to the left. The bottom 9 images show 
gas-temperature maps with similar physical scales from a simulation of a single BCG in which precipitation of cold clouds out of the CGM 
triggers jets from the central AGN Li & Bryan (2014a, b). Cold gas in the simulations (show in blue) displays structures similar in extent 
and axis ratio to the UV structures seen in the star-forming CLASH clusters. Time stamps in the simulation frames show the time in Gyrs 
since the beginning of the simulation. (The first major AGN outburst occurs at ~ 0.3 Gyrs.) 
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Figure 2. Ultraviolet (280 nm rest-frame) broad-band HST images of the central 50 kpc X 50 kpc of 16 CLASH clusters with moderate to 
high entropy cores, smoothed with a gaussian having FWHM of 0.3”. Most of the BCGs in these clusters do not show much 2-d structure 
beyond a faint, central enhancement. Exceptions are the ml206 BCG, with a tadpole-shaped central feature that is likely to be the lensed 
light of a background galaxy, and the ml 149 BCG which has a prominent, lensed spiral galaxy just SE of its core. All of the BCGs are 
centered in each image. In some cases the UV source is very compact, likely to be a nuclear star cluster or a weak AGN. Images are oriented 
with North up and East to the left. 
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Table 3 

UV and NIR Morphological Parameters 


Name 

Aperture 

Centroid Diff b 

dz 

PA (UV) a 

PA (IR) a 

PA Diff b 

dz 

Axis Ratio 

db 

Axis Ratio c 


(arcsec) 

(arcsec) 

(arcsec) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(UV) 


(IR) 

ml931 

4.88 

0.83 

0.01 

-35 

-31.38 

3.2 

0.4 

0.61 

0.01 

0.83 

rl532 

4.88 

0.06 

0.01 

41 

52.91 

12.2 

0.9 

0.79 

0.01 

0.79 

m0329 

4.88 

0.23 

0.02 

-49 

-5.36 

43.3 

1.4 

0.76 

0.01 

0.94 

m0429 

4.88 

1.12 

0.02 

18 

-6.50 

24.7 

1.8 

0.75 

0.01 

0.83 

ml423 

4.88 

0.09 

0.04 

21 

29.07 

7.6 

3.5 

0.65 

0.05 

0.77 

rl347 

4.88 

0.24 

0.03 

-49 

-6.88 

42.2 

10.3 

0.93 

0.02 

0.88 

mlll5 

4.88 

0.33 

0.04 

-43 

-40.29 

2.6 

2.4 

0.78 

0.02 

0.86 

ms2137 

2.28 

0.05 

0.01 

71 

61.08 

10.4 

6.7 

0.96 

0.01 

0.97 

ml720 

4.88 

0.21 

0.04 

-37 

-15.49 

21.8 

34.0 

0.97 

0.06 

0.87 

A383 

4.88 

0.12 

0.05 

27 

13.50 

13.2 

7.2 

0.88 

0.02 

0.96 

r2129 

2.28 

0.19 

0.03 

69 

68.24 

1.1 

3.5 

0.74 

0.02 

0.88 

m0744 

2.28 

0.19 

0.03 

60 

27.52 

32.3 

18.5 

0.91 

0.05 

0.93 

ml311 

2.28 

0.20 

0.04 

-66 

-45.33 

20.8 

7.1 

0.81 

0.05 

0.94 

A2261 

2.28 

0.18 

0.03 

-44 

-27.25 

16.6 

12.1 

0.94 

0.02 

0.98 

ml206 

4.88 

0.13 

0.06 

44 

-77.47 

58.6 

37.0 

1.00 

0.08 

0.76 

A1423 

1.63 

0.15 

0.03 

66 

63.53 

2.9 

6.8 

0.85 

0.03 

0.87 

a209 

1.63 

0.27 

0.03 

-42 

-50.97 

8.7 

6.2 

0.85 

0.03 

0.96 

A611 

1.63 

0.43 

0.06 

38 

36.53 

1.8 

33.5 

0.95 

0.12 

0.95 

cl226 

2.28 

0.58 

0.10 

-79 

-85.89 

6.5 

11.7 

0.31 

0.17 

0.71 

r2248 

1.30 

0.03 

0.01 

64 

50.34 

13.3 

7.2 

0.94 

0.01 

0.95 


Note. — This table includes the 20 X-ray selected CLASH clusters, ordered as in Table 2. 
a Position angles are reported with the convention of the longest axis pointing East of North, in degrees. The rotation solution 
provided in the text computes the angle of the first component (xx) from the x axis in the image. (For our mosaics, the x axis runs 
east-west.) 

b The columns labeled Centroid Diff and PA Diff contain the absolute magnitude of the differences between the IR and UV centroids 
and the IR and UV position angles, respectively, as measured over identical regions on the sky. Note that while the UV and IR 
centroids are all within an arcsecond, only clusters with prominent star formation show significant (> 3cr) differences in the position 
angles. 

c The axis ratio in the IR has negligible statistical uncertainty compared to the ratio in the UV. The ultraviolet emission appears to 
be somewhat more elongated than the IR for star-forming galaxies; the outlier among the quiescent BCGs is in cl226 at z = 0.89, 
which shows an interestingly elongated UV morphology in its core. 


the far infrared or Ha emission (e.g., Dominguez Sanchez 
et al. 2014; Kennicutt 1998; Calzetti et al. 2010). While 
brightest cluster galaxies (BCGs) are either non-existent 
or very rare in these low-redshift multi wavelength stud¬ 
ies, the existing data on BCGs indicate similar ratios of 
unobscured to total star formation (e.g., Donahue et al. 
2011 ). 

For the handful of BCGs with mid-infrared SFR es¬ 
timates based on detections from Spitzer and Herschel 
(Rawle et al. 2012; Hoffer et al. 2012) the correlation 
between MIR SFR and UV SFR is strong (e.g., Hof¬ 
fer et al. 2012). We note that several of the CLASH 
BCGs have nearly starburst-level luminosities (ml931, 
rl531, m0329) as indicated by unabsorbed UV SFRs of 
> 25 Mq yr -1 . ml532 has a mid-IR based SFR ~ 110 
M o yr _1 (0’Dea et al. 2010; Hoffer et al. 2012), which 
is similar to but somewhat higher than the unabsorbed 
UV star formation rate of ~ 80 M 0 yr _1 . Therefore the 
total (summed) star formation rates for CLASH BCGs 
could be 2-3 times higher than our estimates based on 
UV alone. 

All of the BCGs with unobscured UV SFRs > 1 M 0 
yr -1 reside in clusters with low central gas entropy (en¬ 
tropy defined as K = kTrie 2 ^ 3 ). We note that most of 
the central gas entropies reported in Table 2 are from 
the ACCEPT sample (Cavagnolo et al. 2009), with the 
exception of m0416 and m2129, which did not have pub¬ 
lic observations available at the time of the ACCEPT 
project. Those central gas entropies, which are large in 
both systems, are based on the X-ray gas profiles de¬ 
scribed in Donahue et al. (2014). 


These findings for CLASH BCGs at 0.2 < z < 0.9 are 
in accord with the well-known pattern observed among 
lower-redshift clusters. The fact that the thermodynam¬ 
ical state of the X-ray gas, and particularly the central 
cooling time, is strongly correlated with potentially star¬ 
forming optical emission-line gas was recognized as early 
as Hu et al. (1985). Higher resolution X-ray imaging 
spectroscopy from Chandra allowed improved estimates 
of entropy profiles and central cooling times and led to 
the discovery of a well-defined threshold for multiphase 
gas and star formation at K ~ 30 keV cm 2 or t coo i < 1 
Gyr (Cavagnolo et al. 2008, 2009; Rafferty et al. 2008b). 
If cluster cores do not contain gas below these thresholds 
(which are nearly equivalent in clusters of galaxies with 
kT > a few keV), they do not host BCGs with strong 
radio sources or prominent star-formation signatures. 
The fact that the same pattern holds among CLASH 
BCGs (Figure 3) shows that a similar entropy/cooling- 
time threshold has been in place since at least z ~ 0.5. 

4.2. Morphologies 

The CLASH BCGs show a range of UV morphologies, 
and the nature of the morphology is correlated with UV 
luminosity and star-formation rate. BCGs with little to 
no star formation exhibit rather smooth UV morpholo¬ 
gies similar to that of the underlying stellar population 
of evolved stars. In contrast, the UV morphologies of 
BCGs with UV excesses show clumps and filaments, in¬ 
dicating that active BCGs are forming stars, but not in 
disks or monolithic, single super-clusters. The clumps 
are likely to be clusters of hot stars, while the filaments 
could be unresolved strings of recently formed stars, or 







10 


Donahue et al. 


ep 



2 - * 


10 100 ' 
Entropy (keV cm 2 ) 



Figure 3. Rest frame UV (280 nm)-NIR (1 micron) color of the 
central regions of CLASH BCGs plotted against central ICM gas 

entropy. (Gas entropy here is defined to be K = Tn e 2 ^ 3 , where 
T is temperature in keV and n e is electron number per cubic cm.) 
The connection between the thermodynamic state of circumgalac- 
tic gas surrounding the BCG and star formation in the central 
galaxy is clear. The entropy threshold in the CLASH sample is sim¬ 
ilar to the entropy threshold seen in lower redshift clusters (z < 0.2) 
by Cavagnolo et al. (2009) and others. 


possibly emission-line filaments corresponding to inter¬ 
faces between 10 7 K and 10 4 K gas. These structures 
are similar to those seen in the rest-frame UV images 
of BCGs (O’Dea et al. 2010) and in ground-based blue 
images (McNamara 1997). McNamara (1997) identify 
four structural features: point sources, lobes, disks, and 
amorphous sources. As found in O’Dea et al. (2010) and 
in McNamara (1997), most if not all of the star-forming 
systems in our CLASH sample would fall into the amor¬ 
phous category (with multiple clumps). Only a383 and 
ml720 have linear structures that might be very small 
disks. 

Figure 1 shows UV images of 9 of the 11 systems 
with core entropy < 30keVcm 2 (excluding ms2137 and 
r2129), in order of unobscured star-formation rate. All of 
them show significant structure, and the amount of struc¬ 
ture increases with increasing star-formation rate. The 
most extreme star-forming BCGs exhibit the largest UV 
extents and biggest centroid shifts in comparisons among 
UV and near-IR morphological quantities. The other 
star-forming BCGs are somewhat more compact than the 
most UV-luminous BCGs, and are in general fairly well 
centered on the near-IR images. In these less vigorously 
star-forming BCGs, the UV emission is typically less ex¬ 
tended (smaller-scale moments) than the near-IR light, 
even within a limited aperture centered on the BCG. 
The higher entropy systems in Figure 2 are more uni¬ 
form in appearance, with smooth, symmetric, centrally 
concentrated UV emission profiles. UV emission in these 
systems tracks the the evolved stellar population, which 
dominates the near-IR light emitted by both star-forming 
and quiescent BCGs. 

Table 3 reports morphological parameters that quan¬ 
tify these differences between quiescent and star-forming 
BCGs. UV morphological parameters of quiescent BCGs 
are very similar to those measured in the near-IR images, 
while star-forming BCGs exhibit different axis ratios and 


Figure 4. Centroid and position-angle differences between UV 
and near-IR emission for the 20 X-ray selected CLASH BCGs. Blue 
diamonds show star-forming BCGs. Red squares show quiescent 
BCGs. Error bars show 1-sigma uncertainties. Position-angle dif¬ 
ferences in position angles are usually consistent with zero, given 
the uncertainties, except for the star-forming BCGs with highly 
structured UV morphologies. Also, the objects with the largest 
and most significant centroid differences between UV and near-IR 
emission are star-forming BCGs. 

sometimes different centroids in the UV compared to 
the near-IR. Position angles derived from the infrared 
and the UV images are mostly well correlated. Large 
position-angle differences found in a few of the quies¬ 
cent BCGs are not significant because those galaxies are 
nearly round, leading to large uncertainties in the posi¬ 
tion angles. However, 5 of the 6 BCGs with unobscured 
star-formation rates > 10 M 0 yr _1 have UV position an¬ 
gles that differ from the IR position angles by more than 
3(7. 

4.2.1. Morphological Details of Star-Forming CLASH BCGs 

The most extreme SF BCGs in CLASH are in clus¬ 
ters ml931 and rl532. ml931 has multiple clumps, with 
the largest and most luminous clumps in the center. 
There are several clear filaments, aligned preferentially 
north/south. The knot detail is also visible in contin¬ 
uum images, indicating that the stellar density is also 
enhanced in those regions, not just the emission lines. 
The UV structures extend over 20 kpc NW of the BCG, 
and somewhat to the SE, for an overall elongated ap¬ 
pearance. There is a curved, clearly resolved dust lane 
cupping the southeast quadrant of the clumpy nuclear re¬ 
gion, while the brightest three non-nuclear clumps spray 
to the NE, N, and NW. This BCG exhibits GALEX 
emission at even shorter wavelengths (Table 4). R1532 
also has radial filaments extending > 20 kpc in all direc¬ 
tions, but the longest filaments extend to the south, and 
multiple clumps surround the nucleus and extend south¬ 
ward. Among the CLASH BCGs, this one probably ex¬ 
hibits the most similarity to NGC1275. The brightest 
of its knots has been targeted for HST COS photometry 
and STIS long-slit observations later this year (PI: Don¬ 
ahue.) This BCG also exhibits curved, comma-shaped 
dust lanes wrapping the central knots; dust lanes appear 
alongside filaments and knots. 

M0329 exhibits several clumps and a prominent curved 
filament extending north from a flattened structure. It 
lacks FIR observations but is likely a luminous IR source, 
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given that its unobscured star-formation rate is similar 
to those of other BCG LIRGs. It was also a GALEX 
detection, with a 7.5 arcsecond diameter AB mag of 
21.75 =b 0.22 (NUV) and 22.6 ± 0.6 (FUV). 

M0429 shows what may be two distinct, clearly sepa¬ 
rate, clumps interacting in the center of cluster. Several 
other BCGs have a compact, multi-clump morphology. 
Ml 115 has 3-4 clumps oriented like beads on a curved 
string, with a couple of fainter sources a couple arcsec- 
onds off-center. M1423 has a bright but extended central 
source sitting in an elliptical distribution of fainter UV 
light. It is also a luminous UV source in the observer 
frame, with rest-frame far UV detections of 22.0 zb 0.2 
and 23.1 zb 0.6 AB mag (Table 4). Note that for a BCG 
at 2 = 0.545, the NUV and FUV bandpasses of GALEX 
cover mean rest-frame wavelengths of 147 and 98 nm 
respectively. In the F275W image, rest-frame far-UV 
(A ^ 180 nm), which is the bluest image CLASH ob¬ 
tained for this cluster, there is an elongated (0.3-0.5 arc- 
second) nuclear object located precisely with the BCG, 
oriented about 35 degrees W of N. Its IR estimated SFR 
rate of 46 zb 3 M© yr _1 (Rawle et al. 2012) is 3 times the 
unobscured UV rate in Table 2, an IR excess typical of 
obscured starburst galaxies. 

R1347 is a very well-known active BCG in what was 
once the most luminous X-ray cluster known (Schindler 
et al. 1995) and has a cool core with evidence for a merger 
(Mason et al. 2010). It does not show a spectacular 
UV filament system, at least not as spectacular as other 
CLASH BCGs in this category. But it does have a very 
bright SF region just to the ESE of the nucleus, and a 
faint filament on the opposite side of the nucleus. 

4.3. Comparison with Simulations of AGN Feedback 

Cooling and heating of the circumgalactic medium 
(CGM) around BCGs has been a topic of much inter¬ 
est for 30-40 years (Fabian 1994). AGN feedback is in¬ 
creasingly thought to be the key physical process regu¬ 
lating the CGM (McNamara & Nulsen 2007, 2012), and 
the CLASH clusters with low-entropy cores all appear to 
have BCGs with active nuclei. In the CLASH sample, 
10 of the 11 clusters with Kq < 30 keV cm 2 have po¬ 
tential AGN counterparts in the 1.4 GHz NRAO VLA 
Sky Survey (hereafter NVSS Condon et al. 1998) within 
7 arcseconds of the BCG. The only low-ifo cluster lack¬ 
ing an NVSS counterpart (ml931) has a luminous X-ray 
and radio AGN in its core (Ehlert et al. 2011; Donahue 
et al. 2014). (The ml931 BCG radio source is proba¬ 
bly missing from the NVSS because of blending with the 
more luminous Narrow-Angle Tailed (NAT) radio galaxy 
about 30" south of the BCG. See the VLA radio map in 
Ehlert et al. (2011).) In contrast, only 3 of the other 
13 clusters that are members of both the CLASH sam¬ 
ple and the NVSS sky coverage have NVSS radio sources 
with centroids as close to the BCG. One of these clus¬ 
ters is ml206, which hosts a known lensed feature in the 
core of its BCG, so future high-resolution radio obser¬ 
vations and optical spectroscopy are required to distin¬ 
guish the lensed component from intrinsic radio emis¬ 
sion. The other two are cl226, which has three clumps 
in the UV image of its BCG, and A2261, which shows 
evidence for “scouring” by the merging of supermassive 
black holes in the BCG core (Postman et al. 2012a). The 
only CLASH cluster too far south to be in NVSS’s sky 


coverage, r2248, also lacks a counterpart within an arc 
minute of the BCG in the 843 MHz Sydney University 
Molongolo Sky Survey (Bock et al. 1999; Mauch et al. 
2003, known as SUMMS). It also has only upper limits 
for its far-IR or submm flux from Herschel (Rawle et al. 
2012). Also, the two most extreme star-forming BCGs 
in CLASH show evidence in deep Chandra observations 
for X-ray cavities near the central radio source (Ehlert 
et al. 2011; Hlavacek-Larrondo et al. 2013). 

Recent hydro dynamical simulations of precipitation- 
regulated feedback in BCGs Li & Bryan (2014a, b) pro¬ 
duce star-forming structures of cold gas with extents and 
morphologies very similar to the structures we observed 
in UV images of star-forming CLASH BCGs. In these 
simulations, thermal instabilities produce cold clouds 
that rain into the central AGN, which responds by pro¬ 
pelling bi-polar jets into the ambient medium. Interest¬ 
ingly, the circumgalactic medium initially responds to a 
jet outburst by producing additional cold gas, as the jets 
uplift low-entropy gas from the central region of the BCG 
to higher altitudes and cause it to condense into cold, 
star-forming gas clouds. Eventually the cold gas relaxes 
into a disk, but after that disk of gas is consumed by 
star formation, the AGN reignites as the cooling cycle 
re-establishes itself (Li et al. 2015). In these simulations, 
gas flows in and out of the central region of the galaxy, 
and the stars preferentially form along the jet axis with 
a radial extent similar to what we see in the UV images 
of CLASH BCGs. 

Figure 1 shows a montage of simulated temperature 
maps beneath the montage of UV images of CLASH SF 
BCGs. The coldest gas, which shows up as blue in these 
maps, indicates where stars would be forming in the sim¬ 
ulation. We have selected frames from a movie of a single 
simulated cluster, in which the time indicator in the up¬ 
per left corner gives the time since the beginning of the 
simulation. The first major AGN outburst in this simu¬ 
lation happens ^ 300 Myr after the beginning. Frames 
were selected to be similar in appearance to the CLASH 
UV images in the upper part of the figure. Our inten¬ 
tion is to illustrate the various morphologies that can 
emerge during a single episode of AGN feedback and to 
show how similar they can be to the UV morphologies of 
star-forming BCGs. A movie of this simulation can be 
viewed at https: //vimeo. com/84807876. A similar pic¬ 
ture of circulating, heating and cooling gas, moderated 
by the influence of X-ray cavities in the hot atmosphere 
of groups and clusters is discussed in Brighenti et al. 
(2015). 

Interestingly, Lauer et al. (2005) infer similar 
timescales and outflow phenomena from the frequency 
and morphology of dust structures revealed by HST op¬ 
tical imaging of their sample of 77 elliptical galaxies. The 
dust structures seen in imaging at similar wavelengths in 
CLASH BCGs are not identical to the UV structures, 
but show similar extents. We will provide further details 
on the dust morphology and resolved structures in an 
upcoming paper and PhD thesis by Kevin Fogarty. 

5. CONCLUSIONS 

The uniquely wide bandpass coverage and spatial res¬ 
olution of the CLASH HST dataset covering 25 clusters 
with z ~ 0.2 — 0.9 provides a unique view into the UV ap¬ 
pearance of brightest cluster galaxies in some of the most 
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massive clusters in the universe. Contrary to the quieses- 
cent, red-and-dead reputation of BCGs, their UV appear¬ 
ances exhibit a wide diversity of morphologies, similar 
to those found in the small sample of extreme cool-core 
BCGs that have been studied in the rest-frame UV at low 
redshifts. The only CLASH clusters to host BCGs with 
detectable unobscured star formation are those with low- 
entropy X-ray emitting gas in their centers: the “cool- 
core” clusters. The BCGs with large star-formation rates 
are also the only ones with knots, filaments, and other 
asymmetric features in their UV morphologies. One ex¬ 
ception is ml206, with a central feature that is likely to 
be related to a strongly lensed galaxy that is also respon¬ 
sible for narrow, radial filaments in optical images. 

UV emission from quiescent BCGs is less bright but 
still detectable, showing a smooth, centrally symmet¬ 
ric distribution with a typical UV 280 - NIR color of 
~ 5 — 5.5, which is consistent with low-redshift stud¬ 
ies using GALEX and the XMM Optical Monitor. Some 
of these galaxies show moderate UV excesses, which are 
probably due to the small but finite intrinsic scatter in 
the UV-NIR color of the older, evolved population of 
stars. The unobscured star-formation rates we estimate, 
after correcting for this UV component from the elderly 
population, may be underestimates of the total star for¬ 
mation rates by a factor of 2 or more. If this underesti¬ 
mate increases with increasing star-formation rate, as in 
other strongly star-forming galaxies, the star-formation 
rates in the most extreme CLASH BCGs likely exceed 
100 Mq yr -1 . 

The entropy or cooling-time threshold for star for¬ 
mation in the CLASH BCGs is similar to that seen in 
low-redshift (z < 0.2) clusters by Cavagnolo and others. 


While we cannot rule out the possibility of a somewhat 
lower entropy threshold ~ 20 keV cm 2 , it is interesting to 
note that none of the high-entropy clusters in our sam¬ 
ple have BCGs with excess UV light or clumpy UV mor¬ 
phologies, suggesting that even though our X-ray data 
are limited, we are not missing significant low-entropy 
cores in the CLASH clusters. Two potential exceptions 
worth following up in detail are cl226 at z = 0.89, which 
is extremely hot and has a BCG showing some indica¬ 
tion of multiple components along the east-west axis, and 
ml206, in order to sort out the lensed contribution of its 
very interesting central UV/radio source. 

The high resolution UV rest-frame images of star¬ 
forming BCGs show a delightful variety of morphologies 
and structures. The scales and structures of star-forming 
regions seen in these UV images are similar to those of 
cold-gas structures in recent simulations of AGN feed¬ 
back in idealized BCGs, suggesting that while mergers 
and interactions may play a role in the appearance of 
these star-forming regions, AGN feedback is a contribut¬ 
ing if not a dominant driver of the morphology and dis¬ 
tribution of star formation in BCGs. 
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acknowledges financial support from NSF grants AST- 
0908390, AST-1008134, AST-1210890 and NASA grant 
NNX12AH41G, as well as computational resources from 
NASA, NSF XSEDE and Columbia University. 


APPENDIX 

.1. Lensed UV Features Coinciding with the BCG 

Two of the BCGs show UV features that are lensed background galaxies, unrelated to the cluster: ml206 and ml 149. 
We discuss these here because they are very prominent UV features in the HST imagery, but likely have very little to 
do with the BCG itself. 

The BCG in ml206 shows two extremely long narrow filamentary features, one curved and one nearly radial, diving 
into the nuclear region. The center of the BCG has a curled tadpole-type morphology, which is still visible even in the 
bluest UV images available. As a relatively high central gas entropy system of Kq ~ 70 keV cm 2 , it would have been 
surprising if its BCG turned out to host new stars. Two sets of authors have identified these features as likely to be 
strongly lensed background galaxies with uncertain photometric redshifts of around z ~ 1.2 — 1.7 (Zitrin et al. 2012; 
Eichner et al. 2013). 

The BCG in the cluster core of ml 149 is blended with a beautifully lensed, z = 1.491 background spiral galaxy, 
identified in archival data by Zitrin & Broadhurst (2009) and confirmed spectroscopically by Smith et al. (2009). The 
Lyman a emission from this galaxy falls into the rest-frame UV of the BCG. 

. 2 . BCG Notes and Relevant Survey Data 

The BCG of m0717 is non-trivial to identify unambiguously since it is neither dominant nor centrally located. The 
candidate we selected is one of the brightest galaxies spectroscopically identified to be a member of the cluster by 
Ebeling et al. (2014), at an RA (J2000) of 07:17:32.518 and Dec (J2000) of +37:44:34.84, with = 0.5415. This 
galaxy is located away the nominal center of the red galaxies from the CLASH list in Postman et al. (2012b) and 
Medezinski et al. (2013). There are no BCG candidates near this position. The galaxy identified as the BCG for this 
cluster by Rawle et al. (2012) (associated with a Herschel mid-IR upper limit) is likely to be a foreground elliptical 
galaxy. For completeness, we list the central locations of all of the BCGs utilized for UV and near IR photometry in 
Table 4. Most of these coordinates are nearly identical to those in Postman et al. (2012). 

We cross-matched the CLASH BCG locations (in Table 4) with the NVSS (Condon et al. 1998) and SUMSS (Bock 
et al. 1999; Mauch et al. 2003) catalogs of radio sources, with the criterion that a source should be aligned within 7" to 
be considered a candidate counterpart. We estimated the rest-frame 1.4 GHz luminosities (W Hz _1 hf 0 2 ) by assuming 
that L v oc z/ -1 , such that Li. 4 GHz = 47i\D 2 (l + z) 2 Fiagjiz, where D = IZsinr/Tl is the distance measure to the source, 
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Table 4 

GALEX EUV and NVSS Radio Fluxes 


Name 

BCG 

RA(J2000) 

BCG 

DEC(J2000) 

NUV 

AB 

FUV 

AB 

E(B-V) 

E(B-V) 

Aatc/v 

mag 

A FUV 
mag 

NUV C 

AB 

FUV C 

AB 

NVSS 

mJy 

offset 

(") 

Log Li.4GHz 
W Hz -1 

a209 

01:31:52.54 

-13:36:40.4 

21.14 + 0.08 

22.04 + 0.07 

0.0109 

0.11 

0.08 

21.25 

22.12 

<2.5 


<23.5 

A383 

02:48:03.40 

-03:31:44.9 

19.8 + 0.1 

18.72 + 0.12 

0.0315 

0.32 

0.24 

20.12 

18.96 

40.9 

1 

24.6 

m0329 

03:29:41.56 

-02:11:46.4 

20.96 ± 0.08 

21.25 + 0.12 

0.0615 

0.62 

0.48 

21.58 

21.73 

6.9 

6.67 

24.7 

m0416 

04:16:09.13 

-24:04:03.0 

> 22.1 

> 21.4 

0.0413 

0.42 

0.32 

> 22.5 

> 21.7 

<2.5 


<24.1 

m0429 

04:29:36.05 

-02:53:06.1 








138.8 

0.71 

25.9 

m0647 

06:47:50.64 

+70:14:54.0 

> 21.6 

> 22.0 

0.1103 

1.11 

0.85 

> 22.7 

> 22.9 

<2.5 


<24.5 

m0717 

07:17:32.52 

+37:44:34.9 

> 22.1 

> 22.5 

0.0582 

0.59 

0.45 

> 22.7 

> 23.0 

<2.5 


<24.5 

m0744 

07:44:52.82 

+39:27:26.9 

> 22.0 

> 24.3 

0.0578 

0.58 

0.45 

> 22.6 

> 24.7 

<2.5 


<24.7 

A611 

08:00:56.82 

+36:03:23.6 

> 23.0 

> 23.1 

0.0573 

0.58 

0.44 

> 23.6 

> 23.5 

<2.5 


<23.8 

mlll5 

11:15:51.90 

+01:29:55.4 

21.17 + 0.08 

21.6 + 0.1 

0.0386 

0.39 

0.30 

21.56 

21.90 

16.2 

1.47 

24.8 

mll49 

11:49:35.69 

+22:23:54.6 

> 21.5 

> 21.0 

0.0229 

0.23 

0.18 

> 21.7 

> 21.2 

<2.5 b 


<24.5 

A1423 

11:57:17.36 

+33:36:39.6 

> 21.6 

> 21.1 

0.0192 

0.19 

0.15 

> 21.8 

> 21.2 

<2.5 


<23.5 

ml206 

12:06:12.09 

-08:48:03.3 

> 21.5 

> 21.0 

0.065 

0.65 

0.50 

> 22.2 

> 21.5 

160.9 

1.94 

26.1 

cl226 

12:26:58.25 

+33:32:48.6 

23.9 + 0.17 

> 25 

0.0187 

0.19 

0.14 

24.09 

> 25.1 

4.3 

4.37 

25.2 

ml311 

13:11:01.80 

-03:10:39.8 

> 22.4 

> 23 

0.0304 

0.31 

0.24 

> 22.7 

> 23.2 

<2.5 


<24.4 

rl347 

13:47:30.62 

-11:45:09.4 

21.57 + 0.07 

21.82 + 0.08 

0.0615 

0.62 

0.48 

22.19 

22.30 

45.9 

1.09 

25.5 

ml423 

14:23:47.88 

+24:04:42.5 

21.3 + 0.3 

> 21.4 

0.0277 

0.28 

0.21 

21.58 

> 21.6 

8 

2.59 

25.0 

ml532 

15:32:53.78 

+30:20:59.4 








22.8 

0.48 

25.0 

ml720 

17:20:16.78 

+35:36:26.5 

> 21.8 

> 21.2 

0.0374 

0.38 

0.29 

> 22.2 

> 21.5 

18 

1.93 

25.0 

A2261 

17:22:27.18 

+32:07:57.3 








5.3 

6.15 

23.9 

ml931 

19:31:49.65 

-26:34:32.8 

19.49 + 0.11 

19.2 + 0.1 

0.1107 

1.12 

0.86 

20.61 

20.06 

70 


<24.0 

m2129 

21:29:26.11 

-07:41:27.7 

> 22.1 

> 21.4 

0.0759 

0.76 

0.59 

> 22.9 

> 22.0 

<2.5 


<24.5 

r2129 

21:29:39.96 

+00:05:21.2 

21.85 + 0.12 

21.28 + 0.07 

0.0402 

0.41 

0.31 

22.26 

21.59 

25.4 

1.76 

24.6 

ms2137 

21:40:15.17 

-23:39:40.2 

22.09 ± 0.05 

21.89 + 0.03 

0.0511 

0.51 

0.40 

22.60 

22.29 

3.8 

2.52 

24.1 

r2248 

22:48:43.96 

-44:31:51.1 

> 21.90 

> 21.4 

0.0122 

0.12 

0.09 

> 22.0 

> 21.5 

<6 a 


<24.2 


a SUMMS 843 MHz limit. All the others are from the NVSS at 1.4 GHz. 

b The BCG in mll49 may have a nearby radio counterpart, offset about 10 // in the sky. M1149 is a merging cluster, so the radio source may be 
related to merger processes or the BCG; it could also be unrelated. 

c GALEX NUV and FUV total fluxes and upper limits reported in these columns are corrected for Galactic extinction, 
in which r is the comoving radial distance coordinate and 7Z r is the cosmological radius of curvature. The beams 
of both of these surveys are quite large, 45". In one case, ml931, we found that the BCG source was outshined by a 
NAT source in the field, and so the NVSS does not record a radio source near ml931, when in fact VLA observations 
showed a 70 mJy source there. We completed this search to show that all of the BCGs with UV excesses and irregular 
UV morphologies as seen by HST have AGN activity from a relatively powerful central radio source, and we present 
the results in this appendix. 

The GALEX (Morrissey et al. 2007) archive increased in size since the ACCEPT-based work on BCGs by Hoffer et 
al. (2012). Therefore, we also cross-matched the BCG locations with the final (GR6/7) GALEX source list, utilizing 
the GALEXview service 8 from MAST 9 , and report the results in Table 4. The fluxes are estimates of the total 
flux of the source in AB magnitudes in the NUV bandpass (227 nm, 62 nm FWHM) and in the FUV (152 nm, 27 
nm FWHM). The angular resolution of GALEX (about 5") is considerably worse than that of HST, so very little 
morphological information can be gleaned for GALEX for our targets, and as with the radio, there is a danger of 
spurious associations owing to confusion between source and backgrounds/foregrounds. We note that since GALEX 
was an ultraviolet telescope, the central wavelength of its NUV bandpass for CLASH clusters sits almost 100 nm 
farther to the ultraviolet. For about 1/2 of the sample, Lyman-<a is included in the NUV bandpass of GALEX so 
any SFRs based on these numbers would be quite uncertain. (Since the blue continuum of stars is approximately 
flat in L v to 150 nm, a similar AB magnitude for continuum measured at 170 nm would correspond to a similar 
star formation rate measured at 280 nm (Kennicutt 1998), modulo an intrinsic reddening uncertainty.) Nevertheless 
we note that given the very different effective apertures and wavelengths, there is consistency with the upper limits 
and detections by GALEX with the rest-frame UV 280 nmphotometry we have obtained from the HST data. We 
applied approximate Galactic reddening corrections, following our procedure in Hoffer et al. (2012); we report the 
fluxes corrected for Galactic reddening in Table 4. Three-sigma upper limits for exposures up to 500 seconds were 
estimated from the relations estimated Hoffer et al. (2012). Upper limits for BCGs lacking cataloged detections in 
longer GALEX exposures were estimated based on the fluxes of the faintest detected GALEX catalog sources within 
10 arcminutes of the source with flux uncertainties < 0.35 AB mag (~ 3a). 

Facilities: Hubble, Chandra 
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